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We present a systematic first-principles study of the structural and vibrational properties of 
perovskite-structure EuTiOa. Our calculated phonon spectrum of the high symmetry cubic struc- 
tural prototype shows strong M- and R-point instabilities, indicating a tendency to symmetry- 
lowering structural deformations composed of rotations and tilts of the oxygen octahedra. Subse- 
quent explicit study of 14 different octahedral tilt-patterns showed that the I4/mcm, Imma and 
R3c structures, all with antiferrodistortive rotations of the octahedra, have significantly lower total 
energy than the prototype Pm3m structure. We discuss the dynamical stability of these structures, 
and the influence of the antiferrodistortive structural distortions on the vibrational, optical and 
magnetic properties of EuTiOa, in the context of recent unexplained experimental observations. 



PACS numbers: 77.84.-s, 71.15.Mb, 78.30.-j 

I. INTRODUCTION 

Europium titanate, EUT1O3, was the first ternary com- 
pound of divalent Eu 2+ to be identified* 1 -^ It is a G-type 
antifcrromagnet (AFM) with the low magnetic ordering 
temperature of 5.3 K, consistent with the highly local- 
ized 4/ electrons on the Eu 2+ ions^^ Until recently, its 
structural ground state was thought to be the ideal cubic 
perovskite structure, with Pm3m symmetry.— EuTiOa is 
insulating with a high dielectric constant (~ 400) at low 
temperature, indicating quantum paraelectric behavior 
and proximity to ferroelectric instability^ Strong interac- 
tions between the magnetic and dielectric properties have 
been reported^— including evidence for an unusual third- 
order coupling^ It has been suggested that EuTiOs is 
the prototype for studying quantum paraelectric behav- 
ior in magnetic systems iii 

In addition to these experiments on bulk EuTiOs, 
there has been a series of recent studies on EuTiOs thin 
films 1 ^— motivated in part by the subsequently veri- 
fied prediction that epitaxial strain should induce fer- 
roelectricity and multiferroicityi 16 ' 17 Likewise "chemical 
strain", i.e. doping at the A-site with larger divalent 
ions, has been shown to induce ferroelectricity, e.g. in 
alloys with BaTiOs^ The latter is particularly intrigu- 
ing, as multiferroic (Eu,Ba)TiC>3 ceramics fulfil the re- 
quirements for a solid state search for the electric dipolc 
moment of electron^ 

In this work we re- visit bulk EuTiOs, with a partic- 
ular focus on the possibility that the ground state is 
not, in fact, the high symmetry Pm3m perovskite pro- 
toype, but a lower symmetry structure possessing rota- 
tions or tilts of the oxygen octahedra. Until recently, 
all computational studies of EuTiOs had neglected this 
possibility.— We mentioned recently in Ref. OH that 
in the course of calculating the properties of (Eu,Ba)Ti03 
we noticed a calculated antifcrrodistortcd ground state 



for EuTiC>3. Then in Ref. [T?] the influence of such pos- 
sible oxygen octahedra tiltings on the induced polariza- 
tion in strained EuTiOs was discussed. Finally, a recent 
study-22 confirmed experimentally the presence of a struc- 
tural phase transition in EuTiOs using specific heat mea- 
surements, and suggested, in analogy to SrTiOs, that the 
low-temperature phase should have oxygen rotations and 
1 4 1 mem symmetry. 

Our present ab initio study of EuTiOs was motivated 
by a series of experimental papers ; 23 ' 24 in which infrared 
(IR) reflectivity and time-domain terahertz transmission 
spectra of EuTiOs ceramics yielded intriguing results 
that can be better understood if the crystalline symme- 
try is lower then Pm3m. The studie a 23 ' 24 found three 
polar optical phonons between 6 — 600 K, typical for the 
perovskite structure. Analysis of the oscillator strengths 
(or equivalently, the mode-plasma frequencies) and com- 
parison of their values with those of known perovskite 
prototypes such as PbTi03 (with ferroelectrically active 
A- and B-site cations) and BaTiOs (with active B-sitc 
cations) suggested that (i) The lowest-energy transverse- 
optical (TO) phonon TOl corresponds predominantly 
to the Slater mode,— with opposite vibration of the Ti 
cations and the oxygen octahedra (Fig. Q] left); softening 
of this mode dominates the quantum paraelectric behav- 
ior, (ii) The second lowest-energy T02 phonon consists of 
vibrations of the Eu cation against the TiOe octahedra, 
commonly known as the Last mode, (Fig. [1] center), and 
(iii) The highest frequency T04 (the T03 mode is not 
IR active) phonon represents oxygen octahedral bending, 
usually called an Axe mode (Fig. Q] right). 

However, some features observed in the experiment 
are not understood: (i) A significant temperature depen- 
dence of the oscillator strength was found for the TOl 
and T02 modes, indicating changes in the vibrational 
mode eigenvectors with temperature. The changes are 
most apparent in the temperature range 250 — 400 K (see 
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Figure 5 in Ref . H3) , with the oscillator strength for the 
TOl mode decreasing and for the T02 mode increasing 
with increasing temperature, (ii) The width of the peak 
in the imaginary dielectric function e" that is associated 
with the soft TOl mode increases when temperature is 
lowered. This is unusual, because at low temperature 
anharmonic effects that could give rise to the widening 
of the peak in e" usually decrease (see Figure 3(b) in 
Ref. HH). (hi) The mode at around 420 cm -1 in the 
IR reflectivity spectrum is temperature dependent (see 
Figure 2 in Ref. [241 ). However the Eu2Ti 2 07 pyrochlore 
phase, to which this peak was tentatively assigned, does 
not show a structural phase transition. 

In this work we present a detailed density-functional 
theory based study of the structural and vibrational 
properties of antifcrromagnetic EuTi03. Our main find- 
ing is that the ground state is not in fact the high sym- 
metry Pm3m structure (in analogy with other pcrovskitc 
materials, for example KNbOj^), but a lower symmetry 
phase with tilts and rotations of the oxygen octahedra. 
The calculated properties of the lower symmetry phase 
are more consistent with the measured behavior than 
those calculated for the high symmetry Pm3m structure. 
In addition, we obtain a strong magncto-phonon interac- 
tion in our new ground-state structure. The remainder 
of this paper is organized as follows: In Section [TT] we de- 
scribe the technical details of our calculations; Section iHIl 
contains the results of our calculations of the structural 
and dynamical properties of EuTi03, and a new inter- 
pretation of the IR optical properties along with the in- 
fluence of magnetic fields. We conclude in Section [TV] 



II. CALCULATION DETAILS 

We carried out first-principles density-functional cal- 
culations within the spin-polarized generalized gradient 
approximation (GGA)2i. For the electronic structure 
calculations and structural relaxations we used projec- 
tor augmented-wave potentials as implemented in Vi- 
enna Ab initio Simulation Package (VASP)j^ - — We 
considered the following valence-electron configuration: 
5s 2 5p 6 4/ 7 6s 2 for Eu, 3s 2 3p 6 3d 2 4s 2 for Ti, and 2s 2 2p 4 
for oxygen. 

To account for the strong electron correlation effects on 
the /-shells of Eu atoms, we used the DFT+U scheme^ 2 - 
in Dudarev's approach 3 ^ with an on-site Coulomb pa- 
rameter U=5.7 eV and Hund's exchange J_f/=1.0 eV af- 
ter Ref. We have checked our main conclusions (i.e. 
the stability of the antiferrodistorted phases relative to 
the cubic phase) for several values of these parameters 
(U = 4 to 8 eV) as well as without the LDA+U correc- 
tion and obtained the same result. We used a kinetic 
energy cutoff of 500 eV and a6x6x6 (4x4x4) T- 
centered fc-point mesh for the unit cell (supercell) sim- 
ulations. Spin-orbit interaction was not taken into ac- 
count. The experimental lattice parameter of the cubic 
phase of EuTi0 3 , 3.90 A,M£ was used instead of the 



theoretical value of 3.943 A. For structural relaxations 
which involved changes of the unit cell shape we per- 
formed a fixed- volume relaxation to suppress the known 
GGA overestimation of the lattice volume. In the struc- 
tural relaxations we minimized Hellman-Feynman forces 
to 0.5 meV/A. To investigate the structural stability, we 
utilized 2x2x2 supercells containing 40 atoms. This al- 
lowed us to investigate all possible structural distortions 
which can originate from the zone-boundary instabilities 
in a systematic way and to minimize the error in total 
energy calculations when comparing energies of the dif- 
ferent structures. For each symmetry investigated, the 
initial ideal perovskite structure was distorted according 
to the corresponding oxygen octahedral rotation pattern, 
and then the internal atomic positions and lattice pa- 
rameters were relaxed until the convergence criteria were 
reached. 

The calculations of the dynamical properties were per- 
formed using the force-constant methodi 34 i 35 To investi- 
gate possible dynamical instabilities we used a 2 x 2 x 2 
supercell of the primitive 5-atom EuTiOs unit cell; this 
supercell size allows most common oxygen octahedra tilt- 
ing and rotational patterns, as well as the likely commen- 
surate AFM magnetic orders. The Hellman-Feynman 
forces were calculated for displacements of atoms of up 
to 0.04 A. The dynamical matrix for each q-point in 
the Brillouin zone was constructed by Fourier transform- 
ing the force constants calculated at the r-point and 
the Brillouin zone boundaries. Phonon-mode frequen- 
cies and atomic displacement patterns for each g-point 
were obtained as eigenvalues and eigenvectors of the dy- 
namical matrices. Born effective charges for symmetry- 
inequivalent ions in the simulation cel l 36 i 37 were calcu- 
lated with the Berry-phase technique^ 



III. RESULTS AND DISCUSSION 

A. Properties of the high-symmetry P?n3m 
structure 

In the following we show that assuming Pm3m is the 
ground state structure leads to discrepancies with ex- 
isting experiments in the nature of the long-wavclcngth 
phonon modes and in the mode-plasma frequencies. Be- 
cause of the notorious uncertainty in DFT lattice con- 
stants we investigate wheter these discrepancies can be 
removed by varying the unit cell volume. We find that a 
partial agreement (in the mode-plasma frequencies) can 
be reached, but only if the lattice parameter is increased 
to unphysical values at which the lattice is unstable with 
respect to the polar displacements. Furthermore, our 
analysis of the full phonon spectrum of Pm3m structure 
indicates strong antifcrrodistortive instablities. 

First, we calculated the zone-center phonon frequencies 
and eigenvectors in cubic EuTiOs (see Table [TTJ) and ob- 
tained results consistent with other recent calculations ^ 
As we mentioned above, in Ref. [24| it was suggested that 
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Figure 1. (Color online) Schematic diagrams^ of the atomic 
displacements in (a) Slater mode, in which the B-site cation 
moves in the opposite direction from the rigid oxygen octa- 
hedral cage; (b) Last mode, in which the A-site cation moves 
opposite to the rigid B06 cage, and (c) Axe mode, in which 
the apical oxygens move opposite to the in-plane oxygens and 
the A- and B-site cations do not participate; a deformation 
of the oxygen cage results. 



mode- plasma frequencies for the TOl, T02 and T04 
modes, using the following standard definition^ 
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Here (^m) Qa is the mode-plasma frequency of mode m, 
polarized along direction a, eo is the permittivity of vac- 
uum, Vo is the unit cell volume, (Z£) a/3 is the Born effec- 
tive charge tensor of ion k, rrik is the atomic mass of ion 
fc, and {x m )kR i s the eigenvector of mode m. Eigenvectors 
are normalized to unity: 



(2) 



Table I. . Slater, Last, and Axe mode decomposition (in %) 
of the soft polar mode TOl and second energy polar mode 
T02 for two lattice parameters of cubic EuTi03 perovskites. 



Mode uj 
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Axe 


Cell parameter 


a — 


3.90 A: 








TOl 


60 




60 


37 
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T02 


153 




61 


38 
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Cell parameter 
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3.95 A: 








TOl 


96i 




88 


12 




T02 


123 




27 


72 


1 



the TOl mode is mainly Slater type, based on the ob- 
servation that the measured value of the mode-plasma 
frequency for this mode is in the same range as for typi- 
cal ferroelectric perovskites with active B-site sublattices 
such as BaTi03^ In contrast to the suggestion made in 
Ref. HH, however, our analysis of the mode eigenvectors 
shows that the low-energy TOl mode has a rather high 
contribution from Eu displacements, rather than being a 
straightforward Slater mode. In Table Q] we decompose 
our calculated eigendisplacements of atoms in the TOl 
and T02 modes (calculated at the experimental lattice 
constant) into the Slater, Last and Axe modes discussed 
in Ref. |24j. In addition to providing a comparison with 
the earlier interpretation, such a decomposition provides 
information on the relative activities of the A- and B-site 
cations, and correspondingly the role of the two cationic 
sublattices, in the ferroelectric transition. Our analy- 
sis shows that both the TOl and T02 eigenvectors have 
in fact similar contributions from the Slater and Last 
modes. For the TOl mode we find that the Slater- and 
Last-like displacements are in phase, whereas they are 
out of phase for the T02 mode. This leads to a signifi- 
cant difference in the observed, as well as calculated (see 
below), values of the mode-plasma frequencies for both 
modes. 

Next, to compare with the experimental information 
obtained from IR-rcflcctance spectra^ we calculated the 



We used the Born effective charges of 2.63|e| for Eu, 
7A6\e\ for Ti, -2.18|e| for 0± and -5.74|e| for N which 
we calculated at the experimental lattice constant of 3.90 
A in our calculation of the mode-plasma frequencies. 

Interestingly, our calculated mode-plasma frequencies 
(Table [II]) match well with the experimental values in the 
high temperature (300 — 600 K) range, but are markedly 
different from the low temperature experimental values 
(see Figure 5 of Ref. I24T ). Since our DFT calculations 
do not include thermal effects, particularly the strong 
hardening of the soft mode found with increasing tem- 
perature in EUT1O3, we expect them instead to match 
the low temperature values. As we saw earlier, the soft 
mode eigenvectors show significant volume dependence. 
Therefore, to check whether the discrepancy lies in our 
choice of unit cell volume, we investigated the volume de- 
pendence of the key quantities entering Eqn.[TJ We found 
that the Born effective charges are largely insensitive to 
the volume: for values of the lattice constant a between 
3.8 A and 4.0 A the Born effective charges change from 
2.68 to 2.61 |e| for Eu, from 7.46 to 7.49|e| for Ti, from 
-2.26 to -2.13|e| for O x , and from -5.69 to -5.83|e| for M . 
Therefore, the changes in mode-plasma frequencies must 
be caused by the changes in the mode eigenvectors. 

Next, we varied the lattice parameter to see if we could 
obtain good agreement between the calculated and low 
temperature experimental values of the MPFs. In fact 
we found good agreement only by increasing a to 3.95 A. 
From Table U we can see that at a = 3.95 A the TOl 
mode is predominantly of Slater-type, while the T02 
mode has Last character^ However, the Slater mode 
becomes imaginary, indicating a polar structural instabil- 
ity not observed experimentally in EuTiOa under normal 
conditions. Moreover, increasing the lattice parameter to 
3.95 A is clearly not a physically reasonable choice - the 
lattice constant at low temperature (our DFT results cor- 
respond to K) should be smaller than that at room or 
higher temperature because of thermal expansion of the 
lattice. 

Therefore, we conclude that the observed discrepancy 
between the measured low-temperature MPFs and the 
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Figure 2. Phonon spectrum of ideal cubic perovskite EuTiO:?, 
calculated at the experimental volume with lattice parameter 
a = 3.90 A. The imaginary wavenumbers of the low-lying 
modes at the Brillouin zone boundary M- and R- points in- 
dicate structural instabilities. Ferromagnetic ordering of Eu 
magnetic moments was used for computational convenience. 



calculated values is not due to an inappropriate choice 
of the lattice parameter in our calculations. Next we 
explore further the low-tempcraturc structural properties 
of EuTiOa in attempt to explain the discrepancy. 

Our calculated phonon spectrum of EuTiOg in the 
Pm3m phase (Fig. [5]) indeed reveals strong structural 
instabilities at the R (|,|,|) and M (|, |,0) high 
symmetry points, and along the M-R symmetry line. 
These instabilities are seen as a modes with imaginary 
wavenumbers and are similar to those found in calcula- 
tions for Pm3m SrTiO^^ which is known experimentally 
to have the 1 4/ mem ground state with alternating rota- 
tions of the oxygen octahedra around the c axis (Glazer 
notation^ a°a°c~). As in the case of SrTi03, our cal- 
culated eigenvectors indicate that these instabilities are 
non-polar, and arise from the tilting and rotation of the 
oxygen octahedra. The eigenvector for the M-point insta- 
bility shows in-phase rotations of the oxygen octahedra 
around one or more pseudocubic axes, whereas at the 
R point, octahedra rotate with alternating out-of-phase 
sense. 



B. Search for the structural ground state 

Motivated by our finding of rotational instabilities in 
ideal cubic perovskite EuTi03, we next search for the 
structural ground state by comparing the calculated total 
energies of structures containing different combinations 
of rotations of the TiC>6 octahedra. 

As categorized by Howard and Stokes^ condensation 
of combinations of M- and R-point octahedral rotations 
in the perovskite structure can yield any of 15 possi- 
ble symmetry space-groups (including ideal perovskite 
structure) , summarized in Figure [3l The following space 
group symmetries are allowed: PmZm, 7m3, IA/mmm, 
PA/mbm, 14/ mem, Imma, R3c, Immm, P^/nmc, 
Cmcm, Prima, C2/m, C2/c, P1\jm and PI. We cal- 




Figure 3. (Color online) A schematic diagram, indicating 
the high-symmetry cubic perovskite structure and the pos- 
sible 14 subgroup structures, which are the superpositions 
of oxygen octahedra rotations.— Superpositions are indicated 
by Glazer's notations together with the corresponding crys- 
talline symmetries. The structures indicated by rectangles 
relax to ideal cubic perovskite. The structures indicated by 
the rounded rectangles represent the families of tilted struc- 
tures with energies lower than that of the ideal perovskite 
structure. Arrows indicate that a structure with lower sym- 
metry relaxed to a higher-symmetry one. For example P2\/m 
and its high-symmetry neighbour Pnma structure are both 
unstable and relax to Imma. Likewise the unstable PI and 
C2/c structures relax to the R3c phase. 



eulated the total energies of each of these combinations, 
fully relaxing the ionic positions within the constraints of 
the chosen symmetry, and with the unit cell volume set 
to the experimental value. The shape of the unit cell was 
allowed to change according to the symmetry constraints. 

We found that many of the low symmetry tilt patterns 
relaxed to higher symmetry structures, as indicated by 
the arrows in Figure [3] Three structures remained sta- 
ble - those with i?3c, Imma and 74/ mem space groups. 
All three are stabilized by around 25meV per formula 
unit compared with the prototype Pm3m structure. The 
1 4/ mem is lowest in energy (^27 meV per formula unit 
lower than Pmim), Imma is the second most stable 
structure (~26 meV per formula unit lower), then R2>c 
at ^25 meV per formula unit lower. G-type antifcrro- 
magnetic configuration was found to be of lower energy 
than the ferromagnetic one in all tilted phases. Im- 
portantly, the energy differences between the IA/mcm, 
Imma and i?3c structures are very small, within ^2 meV 
per formula unit, which is the range of numerical errors. 
Therefore, we can not predict which rotation scenario is 
the ground state of EuTiOs. All three structures have 
no imaginary phonons through the entire Brillouin zone 
(see Figure [5]) and so are dynamically stable. Note that 
the IA/mcm structure is obtained by condensing out-of- 
phase rotations of oxygen octahedra about only one axis 
(see Figure [5]), Imma results from the same out-of-phasc 
rotations about two axes, and R3c from out-of-phase ro- 
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Table II. Experimental and theoretical values of the phonon wavenumbers and mode-plasma frequencies (both quantities in 
cm -1 ) for the crystalline structures of EuTiC>3 with Pm3m, I4/mcm and RSc space groups, with ferromagnetic (FM) and 
G-type antiferromagnetic (AFM) magnetic ordering on the A-site cations. Mode-plasma frequency values are given in italics. 
T01, T02 and T04 label the optically active modes in the ideal perovsite structure; these modes are split in the tilted 
structures. 



No. 


Mode 


Exp.*, 


Pm3m 




14/ 


mem 






R3c 








10 K ~ 400 K 


AFM FM 


AFM 


FM 


AFM 


FM 


1 


TOl 


82 1550 125 1250 


67 1289 60 1327 


107 


1323 


101 


1394 


80 


1300 


72 


1353 


2 








128 


1607 


118 


1622 


124 


1500 


117 


1540 


3 


T02 


153 400 159 860 


155 925 153 864 


154 


898 


153 


813 


153 


565 


152 


452 


4 








156 


314 


156 


239 


160 


984 


157 


911 


5 








164 




218 




164 




147 




6 








251 


40 


251 


45 


242 


200 


241 


188 


7 








419 


250 


419 


242 


419 


233 


419 


241 


8 


T04 


539 619 540 625 


537 824 536 820 


523 


732 


522 


729 


514 


737 


513 


732 


9 








531 


718 


530 


715 


533 


733 


532 


728 


10 








802 




805 




803 




804 





a Experimental data from Ref. [24| 



where Bq is the bulk modulus, B' its pressure deriva- 
tive, Vo the equilibrium volume, and Eq is the energy 
minimum corresponding to the equilibrium volume. Our 
calculated parameters for this Eqn., obtained by fitting 
of ab initio total energies, are listed in Table [V] The 
equilibrium values for Vo, E$ and bulk modulus Bq for 
tetragonal, orthorhombic and rhombohcdral phases are 
almost the same, however the pressure derivative B' for 
R3c phase is about 1.5 times lower then for the 1 4/ mem 
phase. This indicates that the internal relaxations in 
the rhombohcdral phase are less sensitive to the volume 
changes than in the tetragonal phase. Also, the i?3c 
structure is more rigid than the cubic and the tetrag- 
onal ones. Measurements of the bulk moduli and their 
pressure derivatives at different temperatures would be 
helpful in distinguishing between these lower-symmetry 
phases. 



C. IR optical properties in structures with tilting 



14/mcm 




I h- I »- 

X z 

Figure 4. (Color online) Schematic representation of EuTiOs 
in 74/ 'mem space group structure with oxygen octahedra tilt- 
ing. Left - view along the z direction, where the out-of-phase 
tilting occurs, right - view along the y direction without any 
tilting. The amount of tilting corresponds to the one obtained 
in calculations. Drawings were produced by VESTA visuali- 
sation software.— 



tations about all three axes. Our calculated structural 
parameters for all phases are collected in Tabic HVl 

To check whether our conclusions on the phase stability 
are affected by volume, we next investigated the athcrmal 
equation of state for the three candidate ground-state 
phases, as well as for the undistorted one. We found that 
the E(V) dependence is adequately described by means 
of the first-order Murnaghan equation of stated 



E — E 



BqV ( (Vp/V) 
B' 



B' 



B'-l 



BqVq 

B'o-l 



(3) 



Finally, we attempt to distinguish between the three 
candidate ground state structures by comparing their cal- 
culated phonon and mode-plasma frequencies with the 
measured optical properties. We find that the calculated 
frequencies of the R3c structure show best agreement 
with experiment, suggesting it as a likely experimental 
structure. 

Our results are summarized in Tables [TT| and IIIII as 
well as in Figure [6] In all cases we used the Born ef- 
fective charges calculated for the Pm3m reference struc- 
ture in our evaluation of the mode plasma frequencis, 
and our reported values do not include the LO-TO split- 
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Table III. Phonon wavenumbers and mode-plasma frequen- 
cies (both quantities in cm -1 ) for the orthorhombic Imma 
phase of EuTiOs, with ferromagnetic (FM) and antiferromag- 
netic (AFM) ordering on the A-site cations. Mode-plasma fre- 
quency values are given in italics. The meaning of the TOl, 
T02 and T04 labels is the same as in Table HU 



No. 



Mode 



AFM 



FM 



9 

10 
11 

12 

13 
14 
15 

16 
17 

18 
19 

20 
21 

22 
23 

24 
25 
26 

27 



TOl 



T02 



T04 



36 
40 

98 
104 
106 
110 
110 
131 

154 
155 
159 

164 

237 
249 
283 

413 
414 

416 
418 

424 
436 
491 
499 

516 
523 
537 

801 



1318 



1233 
1596 



201 
945 



134 
122 



224 



730 
736 
725 



34 
49 

85 
104 
105 
107 
112 
119 

152 
155 
156 

166 

235 
247 
282 

412 
415 

418 
422 

424 
434 
486 
498 

513 
524 
535 

802 



1414 
1479 

1609 

656 
816 
107 

134 
76 



167 
302 



716 
743 
707 



ting. Calculations were performed for the G-type AFM 
ordered phase; the coupling between phonon properties 
and magnetism is discussed at the end of this section. 

First we look at the frequency of the soft TOl mode 
(see Tabic [IT]). Our calculated values are 107 cm -1 and 
128 cm -1 for 1 4/ mem (the symmetry lowering splits the 
TOl mode into one doubly degenerate and one singly de- 
generate c polarized mode), 80 cm -1 and 124 cm -1 (i?3c) 
and 98 to 131 cm -1 (Imma). The experimental value is 
82 cm -1 , suggesting that the R3c structure is most likely. 
We note, however, that the soft-mode frequency is very 
sensitive to the volume of the structure and small changes 
in choice of volume (GGA versus experimental for exam- 
ple) can strongly change the predicted values. Regarding 
the mode plasma frequencies of this mode, the experi- 
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Figure 5. Full phonon spectra of EuTiC>3 in 74/mem (upper 
panel), Imma (middle panel) and R3c (lower panel) struc- 
tures. G-type AFM ordering of Eu magnetic moments was 
used in each case. All three structures are dynamically stable. 
In units of the reciprocal lattice vectors, the high symmetry 
points are as follows. IA/mcm: X = (0,^,0), M = (5,5,0), 
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mental value is 1550 cm -1 , and the average calculated 
values for the three structures are 1417 cm -1 (1 4 /mem), 
1433 cm -1 (R3c) and 1382 cm -1 (Imma) all showing 
comparable agreement with the experimental value. Note 
that the relatively large mode splittings that we obtain 
for TOl modes in the lower symmetry structures could 
explain the experimentally observed large peak width in 
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Table IV. Theoretical structural parameters of EuTiOs 
a°a°c~ (I4/mcm), a°b~b~ (Imma) and a~a~a~ (R3c 
gen octahedral rotation patterns. 

1 4/ 'mem 



for 
oxy- 



Wyckoff 






z 


Eu 2(6) 





0.5 


0.25 


Ti 2(c) 











0(1) 2(a) 








0.25 


O 2 4(/i) 


0.21578 


0.71578 





Cell parameters 








a (A) 


5.494 






c (A) 


7.861 








Imma 






Wyckoff 


X 


V 


z 


Eu 4(e) 





0.25 


0.00221 


Ti 4(6) 








0.5 


0(1) 4(e) 





0.25 


0.54663 


0(2) 8(g) 


0.25 


0.02415 


0.25 


Cell parameters 








a (A) 


5.509 






6 (A) 


7.774 






c(A) 


5.541 








RSc 






Wyckoff 


X 


y 


z 


Eu 2(a) 








0.25 


Ti 2(6) 











O 6(e) 


0.53813 





0.25 


Cell parameters 








a (A) 


5.528 






c(A) 


13.448 







Table V. Calculated Murnaghan equation of state parameters 
for the IA/mcm, Imma and -R3c structures. Vo is the theoret- 
ical volume and ao = (Vb) 1//3 is the theoretical pseudo-cubic 
lattice constant. 





Pm3m 


1 4 /mem 


Imma 


R3c 


Vo (A a ) 


61.29 


61.04 


61.11 


61.08 


a (A) 


3.943 


3.937 


3.939 


3.938 


Bo (GPa) 


162.5 


171.1 


169.4 


179.1 


B'o 


10.82 


7.57 


4.98 


4.55 



the imaginary component of the complex dielectric func- 
tion (sec Figure 3b in Ref. [24T ) when two modes with 
lifted degeneracy are fitted as one mode. 

For completeness, we also list in Tables HT1 and Hill our 
calculated frequencies and mode plasma frequencies of 
the T02- and T04-like modes. There are no striking dif- 
ferences that provide persuasive evidence for one struc- 
ture over another. In the IA/mcm structure, the splitting 
of the T02 mode, which is caused by the symmetry low- 
ering, is small (see modes 3 and 4 in Table Hlf . However 
the mode-plasma frequencies of the symmetry-split T02 



2000 



<*-> 1500 - 
1 

$ 1000 - 



— 

1 



500 



TOl 



• 


Experiment 


X 


Pm3m 


o 


14/ mem 


A 


R3c 


V 


Imma 



1 



T02 



v A 

I V I L__L 



T04 _ 




_L 



100 200 300 400 500 600 
Wavenumber (cm ' ) 

Figure 6. (Color online) Calculated and measured^ mode- 
plasma frequencies, plotted as a function of the mode 
wavenumbers. Our calculated values for four different struc- 
tural symmetries are presented. Values originating from the 
TOl, T02 and T04 modes of the ideal perovskite are grouped 
together; the ungrouped symbols correspond to additional 
modes that appear due to the symmetry lowering. 



modes arc significantly different: Mode 3 consists mainly 
of out-of-phase shifts of Eu and Ti atoms, without oxy- 
gen contribution, whereas the c-axis-polarized mode 4 
has non-zero oxygen displacements which reduce the rel- 
ative shifts of the Ti atoms, resulting in a lower dipolc 
for this mode. The average value of the mode-plasma 
frequencies for these three modes is 703 cm -1 , which is 
significantly lower than the corresponding value in the 
Pm3m structure (925 cm -1 ) but still much higher than 
experimental one (400 cm -1 ). The average value of the 
mode-plasma frequency for the T04 mode (modes 8 and 
9) is 727 cm -1 , compared with 820 cm -1 for the reference 
Pm3m structure. In the R3c phase the T02 mode also 
splits by a small amount (7 cm -1 ) with average mode- 
plasma- frequency of 704 cm -1 , and the T04 mode av- 
eraged mode-plasma frequency is 734 cm -1 ; both values 
are similar to those of the 1 4/ mem structure. Like in the 
I4/mcm structure, mode 1 is built by Eu displacements 
against oxygen cage, with small in-phase contribution of 
Ti sublattice. The eigenvectors are also similar to those 
of the PmSm and 1 4/ mem structures. In the Imma 
structure (see Table ITTT|) . the wavenumbers of the T02- 
group are in the range 154-159 cm -1 , and those of the 
T04-group are in the range 516-537 cm -1 . As in the i?3c 
structure, modes in the range 237-283 cm -1 are activated 
and should be observed as a doublet with a splitting of 
50 cm -1 . 

One interesting observation is the presence of an IR- 
active mode with frequency ~419 enr 1 in all three struc- 
tures (mode 7 in Table HI and 18-19 in Table [TTT]). This 
two-fold degenerate mode (it is a doublet in the low- 
symmetry Imma structure) has a relatively low mode- 
plasma frequency of 250 cm" 1 (and, correspondingly, a 
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low oscillator strength). The eigenvector of this mode 
shows that this vibration originates from the out-of-phasc 
displacements of Ti atoms in the neighboring sublattices. 
However, the atomic displacements are not collinear, 
which leads to an uncompensated dipole in the plane, or- 
thogonal to the c-axis. A mode with a similar frequency 
(^430 cm -1 ) is observed at low temperatures in the ex- 
perimental IR reflectivity data (see Figure 2 in Ref. |24[ ) . 
The fact that experimentally this mode is observable only 
at low temperatures gives an additional indication that 
it could be a fingerprint of the antiferrodistortive phase 
transition. 

To make the comparison with the available experi- 
mental data easier, we show in Fig. [6] the calculated 
and experimental mode-plasma frequencies cLS iX func- 
tion of wavenumbcr. Clearly the calculated values for 
the Pm3m structure are strikingly different from the ex- 
perimental values providing further evidence that this is 
unlikely to be the ground state structure. It remains diffi- 
cult to distinguish between the I A/ mem, R3c and Imma 
structures which all show similarly good agreement with 
current experimental data. Our predicted differences in 
the degeneracies and ordering of the components of the 
TOl and T02 modes, could in principle be distinguished 
experimentally using polarized optical spectroscopy on 
single crystals. 

D. Magneto-structural coupling 

Strong magneto-structural coupling and multiferroic 
effects were predicted previously for Pm3>m EuTiOyi= 
here we investigate whether such behavior persists in 
our newly predicted tilted ground-state structures. In- 
deed, we find that constraining the spins to be aligned 
ferromagnetically decreases the soft-mode wavenumber 
by ^7 cm" 1 for all three tilted structures; this is the 
same as previously calculated for the Pm3m structure^ 
(7 cm -1 ), and close to the value estimated from the low- 
temperature IR reflectivity spectra at different magnetic 
fields (3 cm _1 )i2 We obtain a decrease of about 2 cm -1 
for all modes belonging to the T02 group. The frequency 
of the T04 Axe mode is almost insensitive to the mag- 
netic ordering. The modes corresponding to tilting of 
the oxygen octahedra, on the other hand, show strong 
magncto-phonon coupling, since they cause changes in 
the Eu-O-Eu angles which in turn strongly modify the su- 
perexchange pathways between the magnetic ions. These 
modes are listed in Table [TTI as mode 5. We find that ferro- 
magnetic constraint on the Eu spins hardens this mode 



in the 14/mcm structure from 164 cm -1 to 218 cm -1 , 
whereas in the i?3c structure it is softened by 17 cm . 
Since the tilt mode is not polar, however, the dielectric 
function is not affected. 

Finally we mention that Shvartsman et al^ recently 
measured an unusual off-diagonal magnetoelectric cou- 
pling in EuTiC>3, which they rationalized using our pre- 
dictions of tilted ground states. In their experiments, the 
applied electric field breaks inversion symmetry, allow- 
ing a Dzyaloshinskii-Moriya interaction and subsequent 
magnetoelectric response analogous to that in perovskitc 
FeTiC>3, which also shows magnetic A-site ions and oc- 
tahedral tiltings4^ 



IV. CONCLUSIONS 

In summary, our density functional calculations for 
perovskite EuTiC>3 predict that its ground state consists 
of tilts and rotations of the oxygen octahedra, rather 
than the simple cubic perovskitc structure as previously 
believed. By comparing the total energies of all sym- 
metry allowed tilting patterns, we identified three candi- 
date ground states: a°a°c~ (1 4 /mem), a°b~b~ (Imma) 
and a _ a~a~ (i?3c). We compared the calculated phonon 
properties with available experimental data for the three 
candidate structures to identify the most likely ground 
state. Our calculated energy differences are two small to 
allow us to predict which of these phases is the ground 
state; indeed a phase coexistence might be possible, and 
additional experiments are required for further progress. 
In particular, a search for the infra-red-active mode in 
the range 220-290 cm -1 is a signature of i?3c and Imma 
structures, but is not active in 14/mcm phase. Also the 
magnetic-field dependence of the mode at ^164 cm -1 
would distinguish between the 14/mcm and i?3c struc- 
tures. 
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